Theoretical analysis of magnetic Raman scattering in La2Cu04: two-magnon intensity 

with the inclusion of ring exchange 
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We evaluate the Raman light scattering intensity for the square lattice Heisenberg antiferromagnet 
with plaquette ring exchange Ja ■ With the exchange couplings as fixed before from an accurate fit 
to the spin wave dispersion in La2Cu04, leading in particular to Ja ~ 0.24J, we demonstrate in a 
parameter free calculation that the inclusion of the plaquette exchange contribution to the dispersion 
and the magnon-magnon interaction vertex gives a peak position in Big scattering geometry Emax = 
2.71J which is in excellent agreement with the experimental data. Yet, the intrinsic width and the 
lineshape of the two-magnon remain beyond a descriptions in terms of a spin-only Hamiltonian. 

PACS Numbers: 75.40.Gb, 75.10.Jm, 76.60.Es 
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The magnetic properties of La2Cu04 have been the 
subject of many detailed investigations over the last 
decade. Understanding this undoped parent compound 
of high temperature superconducting cuprates is a pre- 
condition for the many theories which describe metal- 
lic cuprates by doping carriers into a layered antiferro- 
magnet. The conventional starting point for undoped 
cuprates is the two-dimensional (2D) spin- 1/2 Heisen- 
berg model with nearest-neighbor (nn) exchange interac- 
tion J 0. Despite the substantial progress on the the- 
ory of the 2D Heisenberg antiferromagnet j2|, some of 
the experimental facts for La2Cu04 have clearly demon- 
strated that a complete description of the magnetic ex- 
citations requires additional physics not contained in the 
2D Heisenberg model with J only. Examples include the 
asymmetric lineshape of the two-magnon Raman inten- 
sity Q or the infrared optical absorption Q. 

The importance of an additional ring (plaquette) ex- 
change coupling for La2Cu04 recently found direct ex- 
perimental support from the observed dispersion of the 
spin-waves along the magnetic Brillouin zone boundary 
A fit of the experimental results j|] using the theo- 
retical spin-wave dispersion, which consistently includes 
quantum renormalization effects has provided accu- 
rate values for the nn exchange integral J = 1720ii' and 
the ring-exchange coupling Ja — 0.24J. The new es- 
timate for J corrects previous values which have been 
used in past years and which were consistently 10 to 15% 
lower. Also the value for Jn must be considered surpris- 
ingly large, but the spin stiffness and the Neel tempera- 
ture calculated with the new parameter set obtained in 
Ref. were in excellent agreement with the experimen- 
tal data and thereby confirmed the deduced exchange 
coupling parameters. The value of the ring exchange Jn 
also agreed with the strong-coupling expansion studies of 
the three-band Hubbard model in the parameter range 
relevant for Cu02 planes 0. 

The Heisenberg model with J only was previously 
found insufficient to describe the experimentally observed 
asymmetry and width of the Raman spectra of undoped 



cuprate compounds The Raman spectrum of the 
Heisenberg model was investigated within the ladder ap- 
proximation for the magnon-magnon scattering vertex 



many years ago 
row peak locatec 



|J9| and gave an almost symmetric nar- 

at -E'niax ~ S.tiJ. It was shown that 

magnon self-energy and vertex corrections to the ladder 
approximation as well as 4-magnon scattering contribu- 
tions are small and negligible. Numerically, the problem 
was investigated by series expansions around the Ising 
limit H^], exact diagonahzations on small clusters 0, 
and QMC calculations [Q. The results of these stud- 
ies led to a peak position and a two-magnon lineshape 
which were very close to the spin-wave results. Fur- 
thermore, the value of the exchange integral J = 1440K 
for La2Cu04 as extracted from the position of the two- 
magnon Raman peak appears too small in comparison to 
early neutron scattering results for the spin-wave spec- 
trum (J = 1650K), see Ref. ||l^ and references therein. 

The asymmetric lineshape of the two-magnon Raman 
intensity Hhas led to proposals that sp in-pho non in- 
teractions ||l^ , [l5| , [l6|] , resonant phenomena [ |l7yi8|] , purely 
fermionic contributions p9|], or cyclic ring exchange 
pO| , pT 22 1 need to be included beyond the nn Heisenberg 
model. A possible importance of ring exchange for Ra- 
man scattering was conjectured also from numerical cal- 
culations which showed that a finite Jn gives rise to 
additional high energy contributions in the Raman inten- 
sity - yet with little spectral weight. 

The new, accurate estimates for the exchange cou- 
plings therefore demand a theoretical reanalysis of two- 
magnon scattering in La2Cu04, which we investigate in 
the present paper with an emphasis on the effects of ring 
exchange. The question that we address is to what extent 
the new values for the exchange parameters - including 
the finite Jn — consistently describe the two-magnon in- 
tensity in La2Cu04 and whether the Heisenberg model 
with ring exchange alone is sufficient to explain the exper- 
imentally observed lineshape of the Raman intensity. To 
answer this question we calculate the magnon-magnon 
interaction vertex in the presence of the ring exchange 
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term and then solve the ladder equation for the Raman 
scattering vertex. 

We start from the Heisenbere model on a square lattice 
with ring-exchange coupling l2^j2^ , ^ 

H = J Si ■ Sm-s + J' Si • Si^s' + J" Si • S^+a" 
+Ju [(S,.S,)(Sfc.SO 

+(s, •s,)(Sfc-Sj)-(s, •Sfc)(s, -so] (1) 

where J, J', and J" are the first (6), second {6') and 
third {6") nearest neighbor in-plane exchanges. We use 
the Dyson-Maleev representation for the spin operators 



S — ajci 



S+ = V2Sbl, S] = -S + b]b, 



ie A , 



(2) 



(3) 



where A and B denote the two sublattices of the an- 
tifcrromagnet; a|,ai, and are Bose operators and 

5 — 1/2. With the Bogoliubov transformation 

Ok = Ukak + fk/3k 
6k = Uk/Sk + V]^al. 
and by introducing the "coherence factors" 

u^.= [{A^, + E^)/{2E^)]^'^ 

«k = -[(Ak-Sk)/(2£;k)]'/' 

we diagonalize the quadratic part of the Hamiltonian and 
obtain the spin- wave spectrum (for details see H) 



E^. = ^Al^Bl (4) 

Ak = 45[J7 - J'l'il -vt)- JnS^hS + 

-AJ"Sy'{l-4'), (5) 
B^ = 4S{Jj~JuS^jf)4, (6) 

with the momentum dependent coefficients 



z^k = (cosfc^ -f cosfcj^)/2; — cos cos ky-, 
fk — (cos 2kx + cos 2ky)/2 . 



(7) 



Note that the quantum renormalization factors {75} in 
(^ take into account the renormalization of the magnon 
spectrum due to quartic terms ^ . The resulting Hamil- 
tonian reads 



^ = E^k(4«k + /3i^/3k) 



(8) 



k "k.,k^ ■ ^ki,<-^k2,i-^k3,mi?k4,ri : ^ki+k2,k3+k4 



where we used the vector notation 



-Rk = (ak,/5l] 



(9) 



and i, j,m,n — 1, 2. The quartic part in Eq. (|8|) is normal 
ordered, since all "Hartree-Fock" renormalizations are al- 
ready absorbed in the quadratic part of the Hamiltonian, 
cf. Ref. 1^. 

We use the effective Loudon-Fleury ||2^ Hamiltonian 
in Big geometry for the coupling of the incoming and 
outgoing photons to the localized nn spins which are in- 
volved in the two-spin flip Raman process 



Hb = Sj • {Sj+s^ - Sj+sJ : 



(10) 



Sx(y) are unit vectors in the directions x, y respectively, 
and A is a coupling constant which includes the elec- 
tric field vectors of the two photons. For the calcula- 
tion of the Raman light scattering intensity in the ladder 
approximation for repeated magnon-magnon scattering 
processes only the vertex V^k'-k'k ~ ^kk^'k'k needed 



],^. The comphcated vertex function V^'^ 



1,2,1,2 



k2;k3k<i 



was 



calculated with the help of computer algebra. After Bo- 
goliubov transformation the ring-exchange term in Eq. 
(1) gives 1057 (!) contributions with different combina- 
tions of a and /3 operators. The general result for the 
is rather involved, however, for equal 



vertex V, 



(4) 



kik2;k2k4 



momenta it simplifies to 



(4) 



-Bk-k'(^k^k' + E]^E]^') — Ak-k'^k^k' — Ra 



kk':k'k 



2E],E],> 



(11) 



where the effects of the ring-exchange coupling are in- 
cluded in the magnon spectrum (^ and in the additional 
vertex contribution 



Rn — 8Jn 



SkBk'(4 



^k-k' 



2Ak'Sk(!^f, + 
2AkAk'(z^Lk' 



''k,k-kO 



□ 
^k.k' 



k-k',k' 



(12) 
) 



2A^B^,{vi 
2E^E^,{-vi_^,+^l^,)\ 



and z/^ ^, = (cos kx cos k'y + cos ky cos k'^)/2. 

The first two terms in the numerator of Eq. ( pT| ) coin- 
cide with the corresponding result for the 2D nn Heisen- 
berg antifcrromagnct [^||^] but with renormalized coeffi- 
cients 74k, B\^ , and the renormalized magnon dispersion 
i?k in (^)- The time-ordered response function of two- 
magnon Raman scattering is given by 



G{u;) 



dte''^'{T[HR{t)HR{0)]) 



where Hji{t) is the Heisenberg representation of Hr^ Eq. 
( |To| ) . The Raman light scattering intensity is then deter- 
mined by — — (I/tt) ImG(tt'). In the ladder approx- 
imation we obtain the equations [^j9| : 
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(13) 



•4°)(k,c^')rkH, 



•Gj'^(ki,c.i)rk,M (14) 



where G'f''^{]i,uj) = — l/(i?kT'^ — *0+) and /k = cosfc^;- 



cos ky . The resuh reads 



io(^) 



L2{UJ) 



, (15) 



Z?(c.) = 1 + ( J - 2Jn52)^^ + ( J + 6J^S')^ 



-{J -2JaS^){J + QJaS^) 



2 Lq{uj)L2{uj) - Ll{uj) 



where 



LM=^ I ^E^G(")(k,.. + c.')G("'(k,c.') 

(16) 



Ef2 Am 1 

k k 

For Jn = the above formulas reduce to the known re- 
sults of Refs. 

For the numerical calculations we use the parameter 
set determined in Ref. ^ from the fit to the spin wave 
dispersion 

J = 151.9 meV, J' = J" = 0.025 J, Jn = 0.24 J. (17) 

For these exchange couplings the renormalization param- 
eters follow as 



7 = 1.158; 7' = 0.909; 7" = 0.852 
7o° = 2.220; 75 = 1-971; jf, = 1.721 



(18) 



Nevertheless we first explore the Jn dependence of the 
two-magnon Raman intensity as shown in Fig. 1. On 
varying Jo the renormalization parameters (Eq) neces- 
sarily have to be recalculated each time. For J' = J" = 
Jn = the result in Fig. 1 coincides with that ob- 
tained in Rcfs. In a first step we turn on J' and 
J" and observe a shift of the two-magnon peak position 
to lower frequencies and a further downward shift when 
also the the ring-exchange coupling is added. The mag- 
nitude of the peak shift due to the ring exchange is quite 
strong. In addition, due to the finite Jn the spectrum 
near the upper edge of the two-magnon continuum de- 
velops a "foot" structure. For comparison, we plot also 
the Raman scattering intensity for a twice larger ring 
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FIG. 1. Evolution of the two-magnon Raman intensity with 
increasing ring-exchange couplings; different parameter sets 
are indicated in the figure. 



exchange Jp = 0.48J where the above effects are even 
more pronounced. The position of the peak estimated 
with Jq = 0.24J is -Emax = 2.71J and with the abso- 
lute values in ( p^ ) is therefore excellent agreement with 
the experimental result |^ -Bmax = 2. 65 J = 3150cm^^. If 
we turn off the vertex contribution Rn in Eq. (11) and 



therefore keep only the influence of the ring exchange in 
the spin-wave dispersion, we obtain Smax = 2.86 J which 
is still much smaller than the peak position without ring 
exchange (E'max = 3.15J). Therefore, the modification 
of the spectrum due to a finite Jn is the major origin for 
the peak shift to lower energies. 

In Fig. 2 we compare the Raman lineshape, obtained 
for the parameter set (p^ with the experimental data on 
La2Cu04 from Ref. 0]. Although the precise peak posi- 
tion is in excellent agreement with the data, its width is 
only slightly infiuenced by ring-exchange coupling. The 
foot structure near the upper edge of the spectrum does 
result from ring exchange and therefore indeed leads to 
an asymmetry of the lineshape, its weight, however, is 
too small to account for the overall linewidth and line- 
shape. Therefore, other processes must contribute which 
are not described by the Heisenberg model with ring ex- 
change alone. It was proposed, that the damping due 
to spin-phonon coupling ||lj,|l^,|l^ is responsible for the 
peak width which in the current calculation is under- 
estimated by roughly a factor 2. In order to incorpo- 
rate extrinsic sources for damping beyond the extended 
Heisenberg Hamiltonian Eq. (1), we have introduced a 
small damping F into the denominators of L„i{lo)^ Eq. 
(|T6|). As a reasonable choice F is assumed constant over 
the entire Brillouin zone, because it is predominantly the 
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FIG. 2. Comparison of the theoretical resuh without 
(r — 0) and with additional damping (F = 0.13 J) to the 
experimental Raman intensity in La2Cu04 in Big geometry 
taken from Ref. [3]. Calculations were performed with the 
parameter set (16). The three experimental curves belong to 
different incoming photon frequencies; the high-energy inten- 
sity background was substracted. All data sets were normal- 
ized to their corresponding intensity maximum. 



zone boundary magnons which determine the shape of 
the two-magnon peak. We plot the result of this calcula- 
tion for r = 0.13J in Fig. 2. One can see that while the 
width of the experimentally observed spectra can be fit- 
ted in such a way, the asymmetry of the spectrum is not 
accounted for. Since the Raman spectra do depend on 
the incoming light frequency, the two-magnon scattering 
is identified as a resonant process and therefore an addi- 
tional source for the strong asymmetry of the lineshape is 
expected to arise from the coupling to the charge degrees 
of freedom 

In summary, we have investigated the effect of ring 
exchange on two-magnon Raman scattering in the 2D 
Heisenberg antiferromagnet. Using the previously fixed 
exchange coupling parameter set we find in a fit-para- 
meter free calculation an excellent agreement with the 
position of the two-magnon peak in La2Cu04. This re- 
confirms the magnitude of the ring-exchange coupling 
Ja = 0.24J in this material. Ring-exchange creates ad- 
ditional high energy spectral weight in the Raman in- 
tensity, but Jn alone is insufficient to explain the overall 
asymmetric lineshape and linewidth. 
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